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HIGHLIGHTS 


►  Nanoscale  3D  Li-ion  battery  electrodes  were  numerically  reconstructed. 

►  Intricate  3D  battery  architectures,  randomness,  and  connectedness  were  detailed. 

►  The  structure-level  electric  transport  was  studied  on  the  reconstructed  structures. 

►  Simulation  data  are  compared  favorably  with  the  experimental  results. 

►  The  pore-level  electrochemical  transport  study  was  attempted. 
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The  goal  of  this  work  is  to  study  electrical  transport  within  multiphase  structures  of  Li-ion  electrodes 
through  numerical  microstructure  reconstruction  and  direct  evaluation.  Li-ion  battery  electrodes  at 
different  compositions  were  numerically  reconstructed  mimicking  the  experimental  fabrication  process. 
Electrode  material  agglomeration  and  interfaces  between  components,  as  well  as  possible  implication  on 
battery  performance  were  explored.  Distribution  of  electrode  phases  (active  material  particles,  binder, 
conducting  phase  carbon),  and  electrical  conductivity  were  studied  upon  the  reconstructed  3D  micro¬ 
structure.  The  electrical  conductivity  of  the  electrodes  was  estimated  from  the  reconstructed  electrode 
structures,  and  compared  favorably  with  those  measured  experimentally  using  four-probe  tests.  It  has 
been  found  that  macro  mass  transport  properties  of  the  numerically  reconstructed  electrodes  depend 
critically  on  nanoscale  physical  characteristics  of  electrode  components  and  their  contact.  This  work 
revealed  the  localized  nanoscale  electrode  complexity  and  transport  processes  in  great  details,  which  are 
currently  difficult  to  investigate  experimentally  or  by  alternative  models. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

High-capacity  Li-ion  batteries  are  among  the  most  viable 
candidates  for  next  generation,  clean,  and  potentially  fossil  fuels 
independent  electric  vehicles.  Extensive  amount  of  research  has 
already  been  devoted  to  improve  the  performance  of  the  Li-ion 
batteries  through  developing  electrode  materials  [1-5], 
reducing  size  [6]  and  optimizing  shape  [3,7,8]  of  the  electrode 
particles,  and  using  different  additives  [9—11].  Despite  recent 
significant  advances,  sluggish  species  transport— both  electrical 
and  ionic— in  a  complex  electrode  architecture  [12,13],  as  one  of 
the  main  factors  limiting  Li-ion  battery  performance,  has  not 
been  resolved.  Fundamental  understanding  and  quantitative 
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investigation  of  the  electrochemical  transport  processes  in  a  pore- 
level  electrode  microstructure,  as  key  to  resolve  the  sluggish 
transport,  still  remain  elusive.  The  lack  of  detailed  understanding 
has  stagnated  the  development  of  the  next  generation  Li-ion 
batteries. 

The  state-of-the-art  experimental  work  has  routinely  relied  on 
trial-and-error  to  obtain  optimum  battery  electrode  design 
parameters  such  as  thickness,  porosity  [11,14],  and  composition 
[14].  With  these  parameters,  performance  of  the  battery  could  be 
widely  varied  as  the  conductivity  of  the  electrode,  and  also  the 
penetration  of  the  electrolyte  inside  the  porous  electrode  structure 
are  altered.  Prior  work  [15—23]  has  confirmed  not  only  that  local¬ 
ized  microstructural  phenomena  are  important,  but  also  that 
interfacial  interactions  [24-28]  are  critical  in  predicting  the  ulti¬ 
mate  performance  and  stability.  Theoretical  work  of  Li-ion  batteries 
has  focused  on  macroscopic  models  that  do  not  address  the  local¬ 
ized  phenomena  at  nanoscale  [12,13,29-32];  therefore  they  are  not 
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capable  of  addressing  the  intricate  electrode  details,  randomness, 
and  the  effects  that  occur  in  3D  battery  architectures.  Furthermore, 
these  localized  phenomena  cannot  be  properly  accounted  for  in  cell 
performance  with  equivalent  circuit  models  [19,33-39];  nor  can 
they  be  assessed  with  sufficient  efficiency  for  simulation  of  large 
domains,  using  atomistic  [40,41]  or  molecular  dynamics  [42-45] 
simulations.  Microscopic  modeling  approach  could  be  proven 
pivotal  since  the  ignorance  of  electrode  structure  details  has  been 
a  major  impediment  to  investigating  the  processes  occurring  at  the 
pore  scale. 

In  this  study,  we  expanded  our  previous  fuel  cell  work  on  the 
process-based  reconstruction  algorithm  and  direct  numerical 
simulation  (DNS)  approach  [46]  to  study  the  Li-ion  battery  elec¬ 
trode  architectures.  LiFePC^  cathode  has  been  reconstructed 
numerically  mimicking  the  experimental  fabrication  process. 
Detailed  microstructure  study  and  conductivity  calculation  have 
been  performed  on  the  reconstructed  structures  and  compared 
with  the  experimental  results.  This  approach  aims  to  provide 
detailed  insight  into  the  localized  nanoscale  electrode  complexity 
and  transport  processes,  which  are  currently  inaccessible  experi¬ 
mentally  or  by  existing  macroscopic  models.  Battery  evaluation  of 
different  electrodes  was  also  performed  to  investigate  their  charge/ 
discharge  behavior. 

2.  Experimental 

Electrodes  containing  Timcal  Super  P  Li  carbon,  poly  vinyli- 
dene  fluoride  (PVDF,  Kynar,  Elf-Atochem)  as  a  binder,  and  active 
material  LiFePCU  were  fabricated.  The  carbon-coated  LiFePCU 
particles  were  synthesized  in  a  solid-state  synthesis  route  using 
two-step  heat  treatment  up  to  700  °C  in  argon  atmosphere  [47]. 
The  active  material,  binder  and  carbon  (corresponding  to  10  wt%, 
20  wt%,  30  wt%,  and  40  wt%)  were  mixed  thoroughly  in 
l-Methyl-2-Pyrrolidinone  (NMP)  to  form  a  viscous  ink.  This  ink 


Fig.  1.  SEM  image  of  the  cross  section  of  a  Li  ion  battery  electrode  (10%  C,  10%  binder, 
80%  LiFeP04).  The  inset  shows  a  high-magnification  SEM  image. 


3.  Simulation  details 

During  numerical  reconstruction,  a  controlled  stochastic  algo¬ 
rithm,  elaborated  in  our  earlier  work  [46],  has  been  adopted  to 
calculate  the  arrangement  of  each  unit  cell,  i.e.,  location  and 
layering  sequence  of  certain  component  within  a  computational 
domain.  An  electrode  segment  of  500  nm  x  250  nm  x  250  nm  was 
reconstructed  using  unit  cells  with  dimension  of 
5  nm  x  5  nm  x  5  nm  (totally  0.25  million  unit  cells).  Each  unit  cell  is 
assigned  to  a  unique  phase  function,  representing  a  unique 

component  in  the  electrode,  e.g.,  binder  (J[i,j,k)  =  2),  carbon 
(j[i,j,k)  =  1),  active  material  ( f[i,j,k)  =  2),  and  pore  (j[i,j,k )  =  0).  As 
shown  in  Fig.  2,  reconstruction  starts  with  nucleation  or  seed,  i.e., 
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was  coated  on  top  of  a  polyether  ether  ketone  (PEEK)  film  using 
a  doctor  blade,  and  then  dried  in  a  vacuum  oven  at  120  °C 
overnight.  Fig.  1  shows  the  cross-sectional  SEM  image  of  battery 
electrodes  with  10  wt%  carbon  additives.  Total  pore  volume  of 
the  electrodes  was  evaluated  by  isopropyl  alcohol  (IPA) 
impregnation.  Electrode  samples  were  immersed  in  IPA  solution 
for  4  h,  and  then  taken  out  for  weight  measurement  after  wiping 
out  the  IPA  at  the  electrode  surface.  The  porosity  of  the  electrode 
was  calculated  according  to  the  weight  gain  due  to  IPA 
impregnation.  The  electrical  conductivity  of  the  electrode  was 
measured  by  the  four-probe  measurement.  Thickness  of  the 
electrodes  was  measured  from  cross-sectional  SEM  images.  Pt  or 
Cu  wires  with  1  cm  distance  were  used.  DC  currents  ranging 
between  1  and  5  mA  were  increased  step-wise  and  the  potential 
measured  at  the  two  inner  probes  was  used  to  calculate  the 
conductivity. 

Battery  performance  was  evaluated  at  room  temperature  in 
a  homemade  pouch  cell,  with  lithium  foil  as  negative  electrode 
and  coated  LiFeP04  on  Al  foil  as  positive  electrode.  All  the 
cells  contained  1.0  m  LiPF6  in  EC— DMC  (1:1  by  volume)  as 
electrolyte  and  were  assembled  in  an  argon-filled  glove 
box.  Charge  and  discharge  characteristics  were  performed 
between  2.5  and  4.0  V  using  a  multi-channel  battery  tester  from 
MTI  Corp. 


randomly  choosing  seeds  within  the  reconstruction  domain  for 
subsequent  growth,  and  then  followed  by  “growing”  materials 
sequentially.  At  each  growth  step,  a  probability  density  function 
(p(i,j,/<)),  which  is  the  number  of  neighboring  “grown”  cells  is 
computed  for  each  empty  cell  in  the  microstructure.  Then  a  prob¬ 
ability  space,  A,  is  constructed  as  Eq.  (1)  shows,  where  cell  (ij,  k)  is 
the  identifier  of  an  empty  cell.  Note  that  the  frequency  of  cell  (i,j,  k) 
included  in  A  is  its  probability  density  function.  Next,  a  stochastic 
process  is  to  choose  one  cell  (B(i,j,/<))  from  space  A,  i.e.,  B(i,j,k)  e  A, 
as  a  new  cell  to  be  grown.  The  above  process  is  repeated  until  the 
desired  amount  of  grown  materials  is  achieved.  For  electrode 
components  that  only  form  on  the  surface  of  other  preformed 
phases  (e.g.,  carbon  coating  on  electrode  materials),  the  probability 
space  will  only  include  the  empty  cells  located  at  the  surface. 

This  reconstruction  approach  provides  an  effective  way  to 
pinpoint  the  localized  electrical  transport  properties  at  pore  level, 
which  can  be  also  transferred  into  a  structure-level  electrical 
conductivity.  According  to  Fick’s  first  law  with  an  applied  current 
density  (/),  the  structure-level  effective  conductivity  is  derived  as, 

0°e  =  J-U/LEe  [2] 

where  A/  is  the  thickness  of  the  reconstructed  domain  (horizontal 
direction  in  Fig.  2),  and  A Ee  is  the  potential  difference.  The  effective 
conductivity  captures  the  collective  information  of  localized 
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Put  seeds  for  agglomerates  Grow  materials  on  seeds  Coat  carbon  on  materials  Mix  with  carbon  and  binder 

Fig.  2.  Illustration  of  the  numerical  reconstruction  process  of  Li-ion  battery  electrode  in  500  nm  x  250  nm  x  250  nm.  The  electrode  contains  7  wt%  C,  7wt%  binder  and  86  wt%  active 
material  with  a  porosity  of  35%. 


transport  and  strongly  depends  on  the  tortuosity  and  connected¬ 
ness  of  the  reconstructed  structures.  Fixed  electrical  potentials 
were  applied  at  the  front  and  back  boundaries,  e.g.,  0  and  IV, 
respectively.  The  electrical  current  under  this  potential  gradient 
was  determined  by  running  a  full  3D  simulation  in  a  commercial 
computational  fluid  dynamic  (CFD)  software  Fluent  6.3.26  and  used 
to  calculate  the  electrical  conductivity  according  to  Eq.  (2).  Thus, 
a  direct  numerical  evaluation  of  material  and  structure  transport 
properties  of  various  designs  is  possible. 

4.  Results  and  discussion 

Li-ion  electrodes  of  several  compositions  (10  wt%  C,  20  wt%  C, 
30  wt%  C,  and  40  wt%  C)  have  been  numerically  reconstructed.  The 
parameters  used  for  the  reconstruction  are  given  in  Table  1.  The 
electrical  conductivity  of  the  Super  P  Li  Carbon  is  estimated  to  be 
around  300  m/s  which  is  close  to  that  of  carbon  nanotube  yarns 
[48].  Fig.  3  shows  the  four  reconstruction  trials  of  the  electrode 
microstructures  with  10  wt%  carbon  additive,  10  wt%  binder,  80  wt% 
of  LiFeP04,  and  a  porosity  of  40%  measured  by  IPA  impregnation. 
Three-dimensional  microstructures  of  the  electrodes  clearly  show 
the  components  such  as  LiFeP04,  carbon  coating  on  the  LiFeP04 
particles,  carbon  additives,  binder,  and  pores.  Careful  comparison 
between  the  numerically  reconstructed  nanostructure  in  Fig.  3  and 
the  SEM  image  in  Fig.  1  suggests  close  similarity.  Both  images 
indicate  a  random  structure  with  irregular-shaped  agglomerates  of 
LiFeP04,  PVDF  binder,  and  carbon  additives.  The  dark  areas  in  the 
SEM  image  are  voids  between  the  agglomerates  and  the  size  is 
ranging  from  tens  to  hundreds  of  nanometers.  The  numerically 
reconstructed  structures  have  similar  agglomerate  size,  pore  size 
and  morphology.  The  reconstructed  electrode  architecture  is  also 
consistent  with  the  observation  in  the  experimental  literature,  for 


Table  1 

List  of  parameters  used  in  the  numerical  reconstruction  of  the  electrodes. 


Parameter 

Value 

Electrical  conductivity  of  Super  P  Li  Carbon 

300  S/m 

Electrical  conductivity  of  active  material 

1  x  10“7  S/m 

LiFeP04  [54] 

Electrical  conductivity  of  carbon  coating 

1  x  10  4  S/m 

on  active  material  [54] 

Percent  of  carbon  in  the  coating  of  the 

8% 

active  material  particles 

Porosity  of  added  Super  P  Li  carbon 

85% 

Density  of  carbon 

2.2  g/cm3 

Density  of  active  material  LiFeP04  [54] 

3.6  g/cm3 

Density  of  binder  (PVDF)  [55] 

1.76  g/cm3 

Number  of  LiFeP04  nucleation  sites 

15 

Number  of  carbon  coating  nucleation  sites 

50 

Number  of  binder  nucleation  sites 

15 

Number  of  added  carbon  nucleation  sites 

25 

example,  Fig.  5  in  reference  [49],  Fig.  8  in  reference  [50],  and  Fig.  5c 
in  reference  [51]. 

Cross-sectional  microstructures  of  the  electrodes  (10%  C,  10% 
binder,  and  80%  LiFeP04),  in  x-z  and  y-z  planes  from  the  four 
different  trials  (Fig.  3),  are  shown  in  Fig.  4.  It  is  clear  that  the  LiFeP04 
particles  have  a  thin  layer  of  carbon  coating  connected  to  the 
network  of  carbon  additives.  The  cross-sectional  microstructures 
also  show  that  there  are  several  isolated  islands  of  carbon  additives 
separated  by  nonconductive  binders.  These  islands  are  not  con¬ 
nected  to  the  main  carbon  network,  therefore  will  not  contribute  to 
the  electrical  conductivity. 

Fig.  5  shows  the  typical  reconstructed  electrode  structures  with 
20  wt%,  30  wt%,  and  40  wt%  carbon  additive.  Structures  from  only 
one  reconstruction  trials  are  shown.  Since  the  reconstructed 
domain  automatically  tracks  the  interfaces  between  different 
components,  e.g.,  particle  and  carbon  additive,  and  their  interfaces, 
the  phase  connectivity,  particle  surfaces,  and  effective  reaction 
interface  within  the  battery  electrode  can  be  directly  accessed  from 
the  reconstructed  microstructure.  Such  a  statistical  analysis  can 
thus  provide  important  information  on  both  transport  and  reaction 
processes  resulted  from  structural  and  compositional  variations, 
e.g.,  weight  fraction  of  carbon  additives. 

The  impact  of  carbon  weight  fraction  is  two-fold.  First ,  LiFeP04 
particle  agglomeration  varies  with  fraction  of  carbon  additives.  At 
a  constant  nucleation  number,  i.e.,  15,  dispersion  of  LiFeP04  parti¬ 
cles  only  relies  on  the  amount  of  active  materials  which  reduce 
with  increasing  carbons  in  a  fixed  reconstruction  volume.  It  is  seen 
from  the  reconstructed  structures  in  Fig.  5  that  higher  fractions  of 
carbon  additives  result  in  smaller  particle  size.  Though  the  primary 
particle  size  of  LiFeP04  is  around  50-80  nm,  its  agglomeration 
leads  to  size  of  200-300  nm  as  shown  in  Fig.  4.  Second ,  LiFeP04 
agglomeration  impacts  interfacial  area  between  different  phases. 
For  example,  the  interface  between  the  active  materials  (with  or 
without  carbon  coating)  and  carbon  additives,  which  is  critical  to 
electron  conduction  within  the  electrode,  depends  on  agglomera¬ 
tion  of  active  materials  and  distribution  of  carbons.  LiFeP04/C 
interface  is  the  most  active  due  to  highly  conductive  carbon  addi¬ 
tives;  LiFeP04/electrolyte  is  the  least  active  due  to  low  conductivity 
of  LiFeP04.  This  is  because  the  electrical  conductivity  of  carbon 
coating  is  2—3  order  of  magnitude  lower  than  carbon  additives  [52]. 
Surface  reaction  in  Li-ion  batteries  relied  on  how  fast  the  electrons 
can  be  supplied  or  removed  from  the  particle  surface. 

Distributions  of  interfacial  area  between  LiFeP04  in  the  typical 
reconstructed  electrodes  at  different  compositions  are  plotted  in 
Fig.  6.  High  carbon  fractions,  i.e.,  30  wt%  and  40  wt%,  create 
noticeable  interface  discontinuity  (zero  interfacial  area  section) 
along  the  thickness;  while  in  comparison,  reducing  carbon  amount 
tends  to  improve  the  interfacial  distribution  along  the  electrode 
thickness.  Moreover,  the  interfacial  area  per  unit  mass  of  LiFeP04 
increases  with  the  carbon  weight  fraction  as  clearly  demonstrated 
in  Fig.  6.  The  averaged  interfacial  area  per  unit  mass  of  LiFeP04  first 
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Fig.  3.  Numerical  reconstruction  of  Li-ion  battery  electrode  containing  10%  C,  10%  binder  and  80%  active  material  with  a  porosity  of  40%:  (a)  Microstructure  of  electrode,  and  (b) 
carbon  distribution. 


increases  from  13.61  m2/g  at  10  wt%  carbon  to  30.1  m2/g  at  20  wt% 
carbon,  and  then  stabilizes  at  around  60  m2/g  when  the  carbon 
weight  fraction  increases  to  30  and  40  wt%.  This  shows  that  inter¬ 
facial  area  and  therefore  macro-homogeneous  electrode  charac¬ 
teristics  rely  on  the  localized  electrode  structure  and  particle 
agglomeration,  which  serves  as  a  design  consideration  for  better 
activity  and  durability.  This  different  surface  coverage  will  have 
a  significant  impact  on  both  reaction  kinetics  and  electrode  elec¬ 
trical  conductivity,  and  therefore  the  discharge  capacity. 

Figs.  3b  and  5b  shows  arrangement  of  carbon  additives  in  the 
simulation  domains  where  a  continuous  electron  transport  path 
along  the  thickness  of  the  electrode  (normalized  x  direction  in  the 
image)  can  be  easily  identified,  although  at  some  points  carbon 
necks  appear.  These  necks  act  as  electron  transport  barriers  and 
will  contribute  to  conductivity  loss.  When  the  carbon  fraction 
increases  to  20%,  30%,  and  40%,  this  connectivity  increases  which 
ensures  the  electrons  supplied  to  the  LiFeP04  particles  during  the 
electrochemical  reaction.  This  can  be  seen  with  the  distribution  of 
carbon  throughout  the  thickness  of  the  electrode  in  Fig.  7.  Results  of 
only  one  typical  reconstruction  trial  have  been  plotted  for  each 
carbon  fraction.  Fig.  7  shows  that  there  is  a  neck  of  carbon 


connection  in  the  electrode  containing  10%  C  where  carbon  volume 
fraction  is  close  to  zero.  Increasing  carbon  fraction  in  the  electrode 
reduces  the  probability  of  neck  formation  and  ensures  a  continuous 
path  of  electron  transport,  and  therefore  probably  better  electrical 
conductivity. 

Fig.  8  shows  the  effective  electrical  conductivity  vs.  carbon 
weight  fraction  averaged  out  of  8  reconstruction  trials.  Electrical 
conductivity  was  calculated  according  to  Eq.  (2).  It  is  believed  that 
a  domain  size  containing  more  than  20-50  electrode  particles  as 
shown  in  Figs.  2—5  (corresponding  to  particle  size  in  the  range  of 
20-100  nm)  is  sufficient  to  achieve  a  realistically  representative 
electrode  macro  properties,  assuming  there  are  no  major  structural 
defects.  The  experimentally  measured  conductivity  value  is  also 
displayed  in  Fig.  8.  The  results  show  a  good  agreement  between  the 
experimental  and  simulated  values.  Conductivity  clearly  showed  an 
increasing  trend  with  carbon  fraction. 

Fig.  9  shows  the  battery  discharge  performance  of  the  electrodes 
with  different  carbon  contents.  The  discharge  potential  is  similar 
for  different  carbon  fractions  at  the  beginning  of  the  two  phase 
transition  region;  however,  it  starts  to  deviate  from  each  other 
along  the  discharging  process.  Though  the  capacities  of  the 
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Fig.  4.  Cross  sections  of  electrodes  containing  10%  C,  10%  binder  and  80%  active  material  with  40%  porosity,  (a)  X—Z  mid-planes,  and  (b)  Y-Z  mid-planes. 


electrodes  seem  low  (the  main  purpose  of  this  study  is  not  to 
optimize  battery  performance),  there  exists  a  clear  trend  that 
increasing  carbon  content  improves  battery  capacity.  From  10  wt 
%  C  to  40  wt%  C,  the  battery  discharge  capacity  increases  from  93  to 


40wt%C 


Fig.  5.  Influence  of  electrode  composition  (amount  of  carbon  additives)  on  electrode 
microstructure:  (a)  Microstructure  of  electrode,  and  (b)  carbon  distribution. 


125  mAh/g  at  0.1  C,  as  a  result  of  the  strong  effect  of  electrode 
conductivity  on  battery  discharge  capacity  [53].  It  is  well-known 
that  both  electronic  and  ionic  conductivities  are  critical  to  the 
cathode  behavior  of  a  Li-ion  battery.  However,  the  IPA  impregna¬ 
tion  measurement  gave  similar  electrode  porosities— around  40%— 
for  the  four  electrodes  with  different  carbon  weight  fractions  in  our 
study.  This  porosity  translates  to  an  almost  identical  effective  ionic 
conductivity  within  the  electrodes,  i.e.,  0.28  S/m,  for  the  four 
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Fig.  6.  Distribution  of  interfacial  area  between  LiFeP04  and  carbon  additives  along  the 
thickness  of  the  reconstructed  Li-ion  electrodes. 
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Fig.  7.  Volume  fractions  of  transport  carbon  throughout  the  thickness  of  the  recon¬ 
structed  electrodes. 


Fig.  9.  Li-ion  battery  discharge  performance  with  electrodes  of  different  compositions 
at  0.1  C. 


electrode  compositions  using  Bruggeman’s  correction.  Therefore,  it 
seems  that  the  observed  increase  of  battery  discharge  capacity 
results  only  from  a  combined  effect  of  improved  electrode 
conductivity  as  well  as  interfacial  contact  between  active  materials 
and  carbon  additives,  as  clearly  demonstrated  in  Figs.  6  and  8, 
respectively. 

5.  Conclusion 

Li-ion  battery  electrodes  of  different  compositions  have  been 
numerically  reconstructed  to  study  the  nanoscale  three- 
dimensional  microstructure  and  distribution  of  the  components 
throughout  the  electrode.  The  electrical  conductivity  of  the  elec¬ 
trodes  was  evaluated  from  the  reconstructed  electrode  structures, 
and  compared  favorably  with  those  measured  experimentally 
using  four-probe  tests.  It  has  been  found  that  macro  mass  transport 
properties  of  the  numerically  reconstructed  electrodes  depend 
critically  on  nanoscale  physical  characteristics  of  electrode 
components  and  their  contact.  The  reconstructed  multiphase 
structures  helps  to  optimize  the  fabricating  parameters  of  the 
electrode  (porosity,  composition,  extent  of  mixing  of  components 
before  coating  etc).  The  approach  developed  in  this  work  thus 
offers  a  direct  numerical  evaluation  of  transport  properties  in 
structures  of  various  designs. 
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